Control of Coulomb blockade in a mesoscopic Josephson junction using single electron 

tunneling 
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We study a circuit where a mesoscopic Josephson junction (JJ) is embedded in an environment 
consisting of a large bias resistor and a normal metal - superconductor tunnel junction (NIS). The 
effective Coulomb blockade of the JJ can be controlled by the tunneling current through the NIS 
junction leading to transistor-like characteristics. We show using phase correlation theory and 
numerical simulations that substantial current gain with low current noise (i„ < 1 fA/\/Hz) and 
noise temperature (< 0.1 K) can be achieved. Good agreement between our numerical simulations 
and experimental results is obtained. 

PACS numbers: 74.78.Na, 85.25.Am, 85.35.Gv 
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Mesoscopic Josephson junctions (JJ) display interest- 
ing phenomena owing to the conjugate nature of phase 
and charge Q. Coulomb effects cause derealization of 
the macroscopic phase variable across the tunnel junc- 
tion, which leads to the formation of energy bands 0- 
As the band width grows rapidly with the band index, a 
mesoscopic JJ makes it possible to construct novel de- 
vices where the operation is based on controlling the 
transitions between energy levels of the junction, thereby 
controlling the effective Coulomb blockade of the device 
HH- The control is made using a small tunnel current of 
single electrons. Thus, the device distinguishes from the 
ordinary Coulomb blockade devices, like single electron 
transistors (SETs), where the current is adjusted by an 
external, capacitively-coupled voltage 0. 

We have investigated the circuit where a mesoscopic 
Josephson junction (JJ) is embedded in an environment 
consisting of a large bias resistor and a normal metal 
- superconductor tunnel junction (NIS). The J J is bi- 
ased in the regime where the system becomes a two-level 
system with two distinct Coulomb blockade strengths. 
The effective Coulomb blockade of the JJ is controlled 
by the tunneling current through the NIS junction, which 
leads to transistor-like characteristics. This device, called 
the Bloch oscillating transistor (BOT) pj, provides a 
low-noise current amplifier whose input impedance level 
makes it an intermediate device between the ultimate low 
temperature amplifiers, the SQUID and the SET 0. 

In this Article, we present experimental results on the 
basic properties of BOTs, and compare them with com- 
puter simulations based on time dependent phase corre- 
lation theory for electron and Cooper pair tunneling. We 
show that this provides a way to model the devices quan- 
titatively. Noise properties of the devices are discussed, 
the conclusion being that ultra low-noise current ampli- 
fiers (current noise i n < 1 fA/vHz referred to input) can 
be built on the basis of controlled JJs. The simulated 
results are shown to be in a good agreement with our 
experimental findings. 

Schematically our device is shown in Fig. 1 (left 




FIG. 1: Schematics of the circuit (left frame) for controlling 
the Coulomb blockade of Cooper pairs. The injector of sin- 
gle electron current, the superconductor-normal junction, is 
denoted by NIS while the Josephson junction is coined by 
JJ. Symbols for different biasing voltages and the resulting 
currents are also given in the figure. Qi denotes the island 
charge that is governed by Eq. (1). See the caption of Table 
I for other parameters. The right frame illustrates the en- 
ergy vs. charge diagram used in the simulations. The distri- 
bution at the Brillouin zone boundary illustrates the charge 
fluctuations due to the resistive environment. Arrows indi- 
cate possible tunneling processes in the BOT: CPT - Cooper 
pair tunneling, ICPT - inelastic Cooper pair tunneling, SET 
- single electron tunneling, Zener - tunneling up to higher en- 
ergy band. Tcp denotes the tunneling rate for both CPT and 
ICPT. 



frame). A Josephson junction (J J) connects a supercon- 
ducting island to the emitter electrode (E). A Normal- 
Insulator-Superconductor (NIS) junction connects the 
normal base electrode (B) to the island. Furthermore, a 
high-impedance, thin film resistor Rc connects the island 
to the collector electrode (C). When the isolation resis- 
tance Rc Rq = h/Ae 2 , charge fluctuations on the 
JJ are small and the dynamics of the junction consists 
purely of Bloch reflections at the Brillouin zone bound- 
ary, interrupted by occasional Zener tunneling between 
the bands 0. 

Fundamentally, the operation of a BOT is understood 
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TABLE I: Parameters for our samples. Tunneling resistances 
(in kfi) of the Josephson and NIS junctions are given by 
i?y J and Rt' S , respectively. Rc denotes the environmen- 
tal impedance of the Josephson junctions. The division of 
capacitance between the two junctions, Cnis and Cjj, are 
estimated on the basis of the measured resistances and geo- 
metrical dimensions. The last column indicates the minimum 
Ey 1 " and maximum E 1 ] 1 ' 1 * values of the Josephson energy in 
/j,eV. 



on the basis of the band structure of a Josephson junc- 
tion. The energy versus (quasi)charge Qi of a J J is shown 
with a thick uniform curve in the right frame of Fig. 
1. When the collector is voltage biased by Vc, the "is- 
land" charge Qi tends to relax through Rc towards the 
value Qi = (Cjj + C NIS )V C - C NIS V B ~ C S V C , where 
Cs = Cjj + Cnis- At Qi = e, Cooper pair (CP) tun- 
neling returns the system back to the state Qi = —e. 
Repetition of this Bloch oscillation cycle produces a net 
current from the collector to the emitter. The competing 
process for the CP tunneling at Qi — e is Zener tunneling 
up to the second energy band, where Bloch oscillations 
do not occur and the current flow stops (Coulomb block- 
ade voltage dE/dQi > Vc)- From the second band the 
system may return back to the first one by single-electron 
tunneling through the base junction. Thus, as a result of 
a single electron tunneling event through the NIS junc- 
tion, a sequence of Cooper pair pulses is injected through 
the J J, which leads to a current gain equaling the aver- 
age number of electrons in the CP sequence. This kind 
of action is below referred to as regular BOT operation. 

The above picture is, however, an overidealization in 
most experimental cases. In a real situation all in- 
put electrons do not produce output current pulses, but 
rather cause intraband transitions. It is also hard to fab- 
ricate large enough thin film resistances to achieve coher- 
ent Bloch oscillations. Qualitatively, the picture of the 
two-level system still works, but the measured current 
gains are clearly less than the calculated values . When 
Rc becomes relatively small, ~ WRq, charge fluctua- 
tions grow and inelastic Cooper pair tunneling may take 
place clearly before the Brillouin zone boundary. This 
makes the dynamics more complicated and our present 
experiments and simulations are aimed at clarifying this 
question. 

In our computational model the island charge as func- 
tion of time is obtained by integrating the equation 
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where the first term represents the charge relaxation 
through the collector resistor, and the last two terms 



represent tunneling current in the NIS and J J, respec- 
tively. {dQi/dt) NIS contains only quasiparticle current 
while (dQi/dt) 7 7 includes both quasiparticle and Cooper 
pair tunneling. The base junction is current biased, as 
in the experiments, but the gate capacitance Cb ~ 1 pF 
see Fig. 1) converts it effectively into voltage bias Vb 



9] 



To integrate Eq. 0) in the presence of finite Josephson 
coupling and electromagnetic environment, we compute 
the Cooper pair tunneling through the J J (Tcp) and 
quasiparticle tunneling through both junctions (T^ J 



urn 



Tqp) using the phase-fluctuation P (E')-theory 
As a modification to the standard theory, we use time 
dependent voltages Vi and Vi — Vb across the JJ and 
NIS junctions, respectively. Vi is given by 



Vi = {C N is/Cv)V b + Qi/C x 



(2) 



A basically similar modelling approach has been em- 
ployed by Kuzmin et al. 12| when investigating the 
role of Zener tunneling on the IV-curves of ultra small 
Josephson junctions. The difference is that we include 
the effect of quantum fluctuations in our model on top of 
the thermal noise which was employed in Ref. |l3| only. 

According to the P (-E)-theory, tunneling does not hap- 
pen strictly at Qi — e, but is rather represented with a fi- 
nite distribution, which is schematically shown in Fig. ^ 
Within this picture, one can still use the concept of band 
structure, the interpretation being that values Qi < e 
correspond to the first band, and Qi > e to the second 
band. The band gaps, where Zener tunneling takes place, 
are now reflected in the probabilities at which the junc- 
tion may pass from the lower to the higher band (see 
below) . 

The most critical assumption is made while comput- 
ing the Cooper pair tunneling rate, for which the lowest 
order theory is valid when EjP (2eV) <C 1 0. For the 
samples analyzed in this paper, max (P (2eV)) « 0.3/ Ec 
and max(-Ej-) w 1.7 Ec, so that second order effects may 
be expected to be small. The Cooper pair tunneling rate 
can then be computed as 



T C p (Vi) 



(3) 



where the function P (E) is defined as 



P(E) = "hh I dt ^(j{t)+^-Et) . (4) 



The phase correlation function J(t) = ([<p(t) — Lp(0)]<p(0)) 
which takes into account the fluctuations of the phase 
ip(t) on the junction. J(t) is calculable from the real part 
of the environmental impedance [lll |. Here we have also 
made the assumption that the effect of the environment is 
exclusively due to the collector resistance for which 
the P(E) function was calculated numerically. 
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In the limit R c / ' Rq > E c /kT, the CP tunneling 
rate as given by Eq. [3] becomes a Gaussian distribu- 
tion centered around Vj — e/Cs [ll|. If furthermore 
Ec/kT ^> 1, the distribution narrows to a delta spike 
T C p = (irE 2 j/2h)6{2eV -4E C ) which equals the equa- 
tion obtained from the band model by neglecting the ef- 
fect of the environment |l5j. Hence, the basic features 
of the band model are embedded in the peaked tunnel- 
ing probabilities. In a proper band model, however, the 
capacitance of the JJ entering Eq. (JJJ would be nonlin- 
ear and given by (d 2 E / dQ 2 )" 1 . In our model, we take 
the capacitance as constant, which is valid in the limit 
Ej — ► only. 

Our sample parameters are given in Table I. Details of 
sample manufacturing and experimental techniques can 
be found in Ref. 0. Measured and computed Ic Vc- 
curves for sample I with various Ej are shown in Fig. 
2, each at several values of Ib- When Zg = 0, a weak 
Coulomb blockade is visible at zero bias. The peak in 
current at non-zero voltages reflects Cooper pair tunnel- 
ing processes, which are enhanced by single-electron cur- 
rent through the base. The current gain j3 is found to be 
maximized in the region with negative slope: (3 max = 3.2 
at V < (I B = 3.3 nA) and (3 max = 3.0 at V > 
(Ib = 2.5 nA). From our simulations we get (3 ma x = 2.8 
(Jb = 1.7 nA) and 2.5 (7s = 1.7 nA), respectively. The 
relatively small maximum gains are caused by large cur- 
rent fluctuations owing to Rc — 23 kil. Fig. 3 displays 
IcVc-cmves for sample II with Ej/Ec = 1.7. Both the 
measured data 0] and the simulated IcVc-cuives dis- 
play hysteretic behavior. The computed curves in the 
regular BOT regime are seen to display growing hystere- 
sis with increasing base current, in accordance with the 
measured data. Especially with large Ej, the current 
peak is at slightly lower voltages in the computed curves 
than in the measured data indicating that the computed 
probability of Zener tunneling is larger than in reality. 
This is a sign that our simplified picture of the energy 
band structure fails for large Ej. 

The nonsymmetrical nature of the curves in Fig. 2 
suggests that the mechanism for current gain is different 
at opposite biasing polarities. This is illustrated in the 
time traces of Qi in Figs. 4a and 4b. In the regular bi- 
asing case (Vc < 0, Fig. 4a) single-electron tunneling is 
seen to drive the system downwards from the upper band 
and, thus, it tries to restore Bloch oscillations, i.e. the 
dynamics is essentially similar to the original BOT oper- 
ation even though incoherent tunneling phenomena mod- 
ifies it strongly. When Vc > 0, on the contrary, Ig tends 
to drive the system into the second band, which leads to 
a suppression of 2e-oscillations as seen in Fig. 4b. The 
recovery of 2e-oscillations now takes place through inco- 
herent 2e-tunneling, the tunneling probability of which 
is determined by the "tail" of P (-E)-function. 

Time domain plots of Qi for the two hysteretic Ic- 
branches (see Fig. 3) are shown in Figs. 4c and 4d. On 
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FIG. 2: Sets of IcVc-curves measured (left) and computed 
(right) at T = 90 mK on sample 1. Josephson coupling has 
been varied with Ej/Ec = 1.2 at the topmost, Ej/Ec = 0.7 
at the middle and Ej/Ec = 0.35 at the lowest set. In each 
set the base currents I B =0, 0.4, 0.8, 1.2, 1.6, 2.0, 2.4, 2.8, 3.2, 
and 3.6 nA (in order from bottom to top). The topmost set 
is offset by -20 nA and middle set by -10 nA for clarity. 
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FIG. 3: ZcVc-curves for sample 2 with Ej/Ec = 1.7 at base 
currents Ib=0, 0.4, 0.8, 1.2, and 2.0 nA. The upper set repre- 
sents data measured at T = 100 mK 6]. The lower set is the 
result of simulation using the experimental parameters. It is 
offset by 10 nA for clarity. 

the upper branch (Fig. 4c) the dynamics is again regu- 
lar BOT dynamics, essentially similar to that of Fig. 4a. 
Now it is obvious that the system almost never relaxes to 
its stationary state, since the inverse of RcC^-time con- 
stant is small compared to the tunneling rate Tnis- Sin- 
gle electron tunneling still clearly enhances Cooper pair 
tunneling. At the lower branch, the system seldom re- 
turns to the lowest band, and the base current Ib mostly 
consists of tunneling events causing intraband transitions 
only. Therefore, Cooper pair tunneling is less likely and 
the collector current remains small. Two stable solu- 
tions can coexist since Ib oc (Vb — Vjj). Now, due to 
the current biased base electrode, Vb and Vjj can each 
dynamically assume two different average values, while 
their difference remains the same leading to equal base 
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FIG. 4: Time traces of charge Qi at the operating points de- 
noted by solid circles in Figs. 2 and 3: a) regular BOT oper- 
ation, b) inverted BOT operation, c) regular BOT operation 
in the hysteretic region, and d) operation in the hysteretic 
region with only intraband transitions induced by Ib- The 
arrows indicate single electron tunneling events. 



currents in both cases. 

In our simulations, the transition from non-hysteretic 
to hysteretic behavior can be crossed by increasing Ej. 
This is because increasing Ej enhances the stability of 
the type of solution presented in Fig. 4(c) by increasing 
the rate Tqp according to Eq. ©. Near the hystere- 
sis point, the ratio of base tunneling currents, Ib 1 /Ib 2 i 
made of contributions causing either interband transi- 
tions (/bi ) or intraband transitions (Jg 2 ) is found to vary 
steeply, for example as a function of Ig ■ Consequently, a 
small change in the input current induces a large change 
in the output current, because of the conversion of base 
current between the types Ib x and Ib 2 ■ Large gain may 
then be obtained with ultra low noise. The output noise 
current is independent of source impedance if Cb C-£. 
In this case the noise temperature can be written 



Tn — Zii 



f] 2 



(5) 



where Z in denotes the input impedance and i nout j (5 is 
the output current noise converted to the input. Near the 
hysteresis point, j3 may grow without limit, which makes 
it possible to reduce Tn — ► 0. Note, however, that 
this takes place at the cost of input impedance ~ (3Rc 
and maximum acceptable input signal amplitude. Ex- 
perimentally, we have found a sample with white noise 
of i n — 10 fA/v^Hz referred to the input. This implies 
a noise temperature T/v smaller by a factor of five com- 
pared to the shot noise approximation i 2 n = 2els- The 
result was measured on a sample having a Josephson cou- 
pling energy of order AEc, i.e. slightly out of the validity 
range of our computational model. In our simulations, 
current noise values of order 1 fA/VHz and noise tem- 



peratures below 0.1K have been reached using experi- 
mentally realizable parameters. 

As mentioned in the beginning, our analysis is based 
on perturbation theory in Ej/Ec and the results are 
strictly valid in the limit Ej — ► 0. Nevertheless, we be- 
lieve that our simulations are principally valid at values 
of Ej ~ 1 . More serious problems occur when the width 
of the lowest band becomes exponentially narrow as for 
samples with Ej/Ec ^4 6]. In this case, the dynamics 
is dominated by transitions between higher bands, rather 
than the two lowest ones. As long as there are two major 
bands involved allowing to use a two band approxima- 
tion, our simulated results will be qualitatively correct. 

In summary, we have shown that novel devices can 
be constructed using interlevel transitions, driven by sin- 
gle electron tunneling. The single electron current can 
be used either to drive the JJ in to the blockade state, 
or out from the blockade. Both methods are seen to 
yield substantial current gain, though at slightly different 
values of bias voltage. Comparison of our experimental 
and theoretical results show that time-dependent P(E)- 
theory can be employed quite successfully to model the 
behavior of such devices. The essential features of BOT 
dynamics were found to be recovered even in the pres- 
ence of relatively strong incoherent Cooper pair tunnel- 
ing and the JJ could still be understood as a two-level 
system, whose switching is controlled by single charge 
tunneling. Furthermore, our simulations indicate that in- 
traband transitions play a significant role in these three 
terminal devices. The devices may give substantial gain 
by conversion of base current between interlevel and in- 
traband types. This leads to noise powers that are sub- 
stantially less than the values obtained from the input 
shot noise approximation. 
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